Acoustic vibrations of anisotropic nanoparticles 
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Acoustic vibrations of nanoparticles made of materials with anisotropic elasticity and nanoparti- 
cles with non-spherical shapes are theoretically investigated using a homogeneous continuum model. 
Cubic, hexagonal and tetragonal symmetries of the elasticity are discussed, as are spheroidal, cuboc- 
tahedral and truncated cuboctahedral shapes. Tools are described to classify the different vibrations 
and for example help identify the modes having a significant low-frequency Raman scattering cross- 
section. Continuous evolutions of the modes starting from those of an isotropic sphere coupled with 
the determination of the irreducible representation of the branches permit some qualitative state- 
ments to be made about the nature of various modes. For spherical nanoparticles, a more accurate 
picture is obtained through projections onto the vibrations of an isotropic sphere. 



PACS numbers: 

I. INTRODUCTION 

The lowest frequency vibrations of isolated nanopar- 
ticles are in the THz range, on the order of the speed 
of sound divided by the dimension. These are commonly 
refered to as confined acoustic phonons and are unrelated 
to optical phonons. There have been numerous experi- 
mental and theoretical studies on the acoustic vibrations 
of nanoparticles in the last few decades. These vibra- 
tions have been observed by a variety of experimental 
techniques including low frequency Raman scattering,A 
time resolved femtosecond pump-probe experiments;^!^ 
infrared absorptiouf^'^ inelastic neutron scattering^ and 
persistent spectral hole burning!^ 

Reasonable estimates for the mode frequencies are ob- 
tained using the 1882 Lamb solution of the continuum 
elastic problem for an elastically isotropic, homogeneous, 
free sphere.® This provided sufficiently good agreement 
to confirm that confined acoustic phonon modes were re- 
ally being observed. However, this model was unable to 
deal with the anisotropy of actual samples. 

Recent advances have permitted the creation of high 
quality elastically anisotropic samples The essential fea- 
tures are (1) a narrow size distribution (2) good crys- 
tallinity so that a significant amount of nanoparticles 
in the sample are mono-domain (3) controlled shape of 
nanoparticles and (4) separation of nanoparticles so that 
they vibrate as independent units. As a result, the vi- 
brational modes of elastically anisotropic nanoparticles 
have been observed. This has created the need for an 
alternative to the Lamb model capable of dealing with 
nanoparticles with lower symmetry. 

In this work, we use the method of Visscher et alJ^ 
which is a standard numerical approach suitable for the 
calculation of the frequencies and the wavefunctions of 
the vibrations of such nanoparticles. The symmetry 
of these modes, their volume variation and their Lamb 
mode parentage are determined and applied to the pre- 



diction of their observation by different experimental 
techniques such as inelastic light scattering. 



II. METHODS 

The situation for an isotropic nanoparticle will now be 
summarized. In this case, the system is spherically sym- 
metric. Thus, vibrational modes can be classified by their 
angular momentum number £ > and its z-component 
m. Modes can also be classified either as torsional (T) 
or spheroidal (S). Finally, modes are also indexed in or- 
der of frequency by n > 1 (n = 1 corresponds to the 
first harmonic (fundamental mode), n = 2 to the second 
harmonic and so on). In the following, we will indicate 
Lamb modes using the compact notation X"^ where X=S 
or X=T. 

All modes can be observed by inelastic neutron scat- 
tering in the typical situation where the wavelength of 
the neutrons is much smaller than the nanoparticle size. 
For nanoparticles whose dimension is small compared to 
the wavelength of light (dipolar approximation) , Raman 
only detects Sq and S2, infrared absorption only detects 
Si"ii and time- resolved femtosecond pump-probe experi- 
ments typically only detect Sq. In the following, we will 
assume that the nanoparticles are small enough so that 
the dipolar approximation holds. 

Nanoparticles with either isotropic or anisotropic elas- 
ticity will be considered in the following. Isotropic elas- 
ticity is considered mainly for comparison with previ- 
ous studies. Anisotropic elasticity is used for perfect 
nanocrystals consisting of a single domain. We will re- 
fer to such nanoparticles as being "mono-domain" in 
the following. In a small nanoparticle, a mono-domain 
structure is not necessarily energetically favorable. For 
example, it is well known that multiply-twinned silver 
nanoparticles are much more stable for certain ranges of 
sizejii 
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A. Calculation of frequencies and associated 
displacements 

The frequencies and their associated displacements for 
an anisotropic nanoparticle have been calculated using 
the approach introduced by Visscher et alr^ which also 
assumes continuum elasticity. Other authors have al- 
ready confirmed that the convergence of this method is 
faster than the convergence of finite element methods at 
least in some cases.-iS, The relevance of continuum elastic- 
ity for nanoparticles has been confirmed using atomistic 
calculationsi^ii^ii^iii for nanoparticles larger than 2-3 nm 
and even for ZnO nanoparticles for which surface relax- 
ation and stress are significant Most of the results pre- 
sented in this paper have been obtained for nanospheres 
whose diameter is 10 nm which is well above these limits. 
It is possible to extrapolate them to different sizes since 
the frequencies vary as the inverse diameter. However 
care should be taken not to consider very small nanopar- 
ticles for which surface effects could significantly alter 
the validity of the continuum approximation. The calcu- 
lational method for the modes gives each mode in terms 
of power series coefficients Uijpqr so that: 

Ui{x,y,z)^^ ^ aijpqrxPy'^z''i (1) 

pqr j=x,y,z 

The power expansion covered 0<p + g + r<20 for 
good convergence for all the modes we are interested in. 
The frequencies for the isotropic spherical case were re- 
produced with very good accuracy. The convergence for 
strongly anisotropic systems is harder to check. Despite 
checking that the frequencies do not significantly change 
when adding more terms to the power expansion, we also 
compared the calculated frequencies with the Finite El- 
ement Mesh Sequence method introduced in a previous 
work.l- 

In this work, all the displacements have been normal- 
ized according to equation [2] where V is the volume of 
the nanoparticlci^ 




B. Group theory 

1. Degeneracy lifting 

In the absence of spherical symmetry, modes are clas- 
sified according to their remaining symmetry. For ex- 
ample, for a spherical nanoparticle with cubic elasticity, 
such as Ag, Au and Si, the system is symmetric under the 
symmetry operations of a cube, which is the 48-clcment 
group Oh- 

Our calculations return a large number of modes which 
have to be considered to interpret inelastic light scatter- 
ing spectra or other experimental results. It is important 



to use not only the frequencies but also the wavefunctions 
in order to do that. Group theory is a very valuable tool 
in this context as it allows for example to identify the Ra- 
man active modes and therefore simplify the assignment 
process. In this work, we will only consider nanoparticles 
whose dimensions are small compared to the wavelength 
of light (dipolar approximation) in order to discuss the 
selection rules of Raman scattering. Table U shows how 
the degeneracy of the Raman and infrared active modes 
of an isotropic sphere- - is lifted or not when lowering the 
symmetry. Only point groups relevant for the rest of this 
paper are considered. 



2. Numerical determination of the irreducible 
representations 

In order to take full advantage of group theory, it is 
important to label the different modes with the corre- 
sponding irreducible representation. Sophisticated and 
specific approaches could be considered to restrict the 
calculations to modes having a well-defined symmetry. 
However, we preferred to keep the numerical approach 
detailed previously because it is more general. We added 
a few calculation steps to determine the irreducible rep- 
resentation from the wavefunctions. It turns out this can 
be achieved very reliably and without much additional 
calculation time. 

S being a symmetry operation of the point group of 
concern, and Ui with i = 1, . . . n a full set of eigenmodes 
having the same frequency, the character of S for this 
irreducible group of vibrations is: 



^{S)=Y, jjj S (u, (r) ) • (S (r) ) <fR (3) 

Such integrals can be calculated accurately and quickly 
when each Cartesian component of the displacement field 
is a sum of terms of the form x^y'^z^ . By calculating 
a few well-chosen characters, it is then straightforward 
to distinguish all the different irreducible representations 
using the character table of the point group. The only 
restriction to apply this procedure is that the degeneracy 
must be known and therefore the convergence must be 
good. Of course, a more general approach is require to 
handle accidental degeneracies. 

In the following we detail some additional ways to 
improve our knowledge of these vibrations. These are 
needed since regarding Raman scattering many modes 
are labeled as Raman active due to the irreducible repre- 
sentation they belong to. However there is not necessar- 
ily an efficient coupling mechanism enabling a significant 
Raman intensity. The following tools are designed to 
somewhat address this problem. 
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Aig 
Aig 
Aig 
Aig 


Aau + Eu 
Aau + Eiu 
Aiu + Eiu 


Eg + Tag 
Aig + Big + Bag + Eg 
Aig + Eig + Eag 
Aig + Eig + Eag 


Aau + Tiu + Tau 
Aau + Blu + Bau + 2 Eu 
Aau + Blu + Bau + Eiu + Eau 
Aiu + Eiu + Eau + Esu 


Aig + Eg + Tig + Tag 
2 Aig + Aag + Big + Bag + 2 Eg 
Aig + Big + Bag + Eig + 2 Eag 
Aig + Eig + Eag + Esg + E4g 



TABLE I: Degeneracy lifting of spheroidal Lamb modes with £ < 4 for various symmetries. The degeneracy lifting for torsional 
modes (with ^ > 0) is obtained by changing the parities {u ^ g) 



C. Volume variation 

The volume of the nanoparticle does not change for ev- 
ery possible vibration. For isotropic spherical nanopar- 
ticles, the volume changes only for the spheroidal £ ~ 
vibration. This volume change is also involved in the 
time-resolved femtosecond pump-probe measurements.'^ 
As a result, it is interesting to calculate it for all the 
vibrations. 

The volume variation corresponds to the flux of the 
displacement through the surface of the particle. The ex- 
pression for a dimensionless volume variation SV is given 
in equation [D Using the divergence theorem, this 2D 
integral can be turned into a volume integral involving 
derivatives of x^y'^z^ functions which arc already calcu- 
lated in the frame of Visscher's method.^" Therefore this 
quantity can also be accurately and efficiently calculated. 
It can be shown that the volume variation is different 
from zero for fully symmetric vibrations only, i.e. for 
Aig vibrations for the symmetries considered here. 



jj U, (i?) ■'^^ ^ jjj {^i 1 1 



(4) 

Of course, the value of 6V depends on the normaliza- 
tion (see equation [2|) . 



D. Smooth variation of parameters 

In order to follow how the vibrations evolve when low- 
ering the symmetry, it is possible to follow the frequencies 
of the different modes while slowly varying the parame- 
ters (for example the shape or the elastic constants of the 
material the nanoparticle is made of). The curve repre- 
senting the variation of one frequency is called a branch 
in the following. Group theory can help plotting such 
branches more reliably because all the different points 
on a given branch share the same irreducible represen- 
tation. Moreover branches having different irreducible 
representations can cross but not branches having the 
same irreducible representation. This originates from the 
coupling between the different branches being due to the 
anisotropy itself whose irreducible representation is Aig 
for the point groups we focus on. As a result, all n*-^ 
modes having a given irreducible representation belong 
to the same branch and can be safely connected. 

In this work, branches due to varying elastic anisotropy 



are calculated by changing the stiffness tensor using 
C{x) = (l-2:)C'"°-ha;C^'" with < a; < 1, C'''° being the 
tensor of the isotropic material (which can be obtained 
by averaging the sound velocities^ or other methods^ 
or using measured longitudinal and transverse sound ve- 
locities) and C™' the anisotropic tensor. The parameters 
used in this work are given in Ref. [2^ . It should be noted 
that such branches are made of fictive materials except 
for x = 1 where the real elastic parameters of the bulk 
material are used. 



E. Projections 

When the lowering of the symmetry is not due to a 
change in the shape of the nanoparticle, it is possible to 
compare the wavefunctions (omitting their time depen- 
dences) of the two systems. As commonly done for atom- 
istic calculationsi^ii^ii^, we calculated the projection of 
the displacements of spherical nanoparticles onto Lamb 

modes. For a given mode whose displacement is Ui ( R 



the projection onto the Lamb mode is defined as: 
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[r)-xI,(^r)£R (5) 



The orthonormality and completeness of the Lamb 
modes implies that, for any i. 



e.m) 



= 1 



(6) 



Likewise, the orthonormality and completeness of the 
modes of the anisotropic nanoparticle implies that, for 
any m and n, 



(7) 



However, a more relevant quantity is obtained by sum- 
ming the squared projections over all the degenerate 
Lamb modes i.e. over m. As a result, the total squared 
projection of Ui onto the subspace spanned by the 2£+l 
modes X? is: 



(8) 
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V (GHz) 



103.3 
106.8 
151.1 
153.0 
161.0 
169.4 
196.4 
196.9 
211.2 
221.8 
232.1 
234.5 
248.9 

331.2 



Eg 
Tau 
Aag 
Tiu 

T2g 
T2g 

Tig 
Eg 
Tau 

T2u 

Aig 
Aig 



Squared Lamb projections 
0.995 Sa 



0.936 
0.996 Ti 
0.499 S3 
0.974 
0.821 Sa 
0.758 T\ 
0.411 
0.838 Sa 
0.748 Ti 
0.466 S3 
0.665 S3 
0.959 si 



+ 0.002 
+ 0.054 
+ 0.002 
+ 0.480 
+ 0.018 
+ 0.160 
+ 0.178 
+ 0.369 
+ 0.113 
+ 0.199 
+ 0.466 
+ 0.127 
+ 0.028 



Si 
Si 

sj 

s; 

% 

S^ + 
Si + 
Si 
SJ 

s; 

Ti 
Sj 



0.930 S?) + 0.039 Si + 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.4 

4.3 



TABLE II: Characteristics of the modes of a mono-domain 
Ag sphere with R = 5 nm are shown, i is the mode index. 
The two largest projections onto Lamb modes are shown. 
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The Lamb projections given in the fourth columns of 
tables im through I VIII correspond to this last quantity. It 
represents the ratio of the energy contained in the pro- 
jection onto Lamb mode X" to that of mode i assuming 
the same frequencies. 



III. APPLICATIONS 



FIG. 1; (Color online) Evolution of the frequency of Raman 
active modes of a mono-domain silver sphere (radius 5 nm) 
with varying elastic anisotropy up to the frequency range of 
the breathing mode. Raman inactive modes are plotted with 
black lines, Eg branches with lines and crosses (red online), 
Tag branches with lines and triangles (red online) and Aig 
branches with lines and full circles (blue online). 



A. Spherical nanocrystals with cubic crystallinity 

1. Spherical mono- domain silver and gold nanoparticles 

Let US first consider the case of a spherical nanopar- 
ticle made of silver. This nanoparticle is mono-domain 
and therefore the stiffness tensor is the same everywhere 
inside the nanoparticle and identical to that of bulk 
silver Table HIl presents the calculated frequencies, irre- 
ducible representations, main projections onto the modes 
of an isotropic silver sphere and volume variations of the 
lowest frequency modes. The six modes i — 1 ... 6 have 
zero frequency and correspond to the rigid rotations and 
translations of the nanoparticle. The branches corre- 
sponding to the lowering of the symmetry when going 
from the isotropic to the anisotropic case are plotted in 
figure [TJ 

Figure [T] clearly shows that the introduction of elas- 
tic anisotropy significantly lifts the degeneracy of most 
modes. One notable exception is the breathing mode 
which corresponds to the second Aig branch. This mode 
is non-degenerate and its frequency hardly changes with 
anisotropy. The other important exception are the dipo- 
lar modes Si which are infrared active and transform into 
Tiu with the same degeneracy. The lowest Raman active 
mode S2 which has degeneracy 5 is split into the low- 
est Eg and T2g branches {i = 7,8 and i = 18, 19, 20 in 



table |TT] respectively) as expected from table HI These 
anisotropic modes have a dominant projection onto 83 
confirming their Lamb mode parentage. 

Similar calculations have been performed for mono- 
domain gold nanoparticles using the elastic constants 
from Ref.[2J and the results are presented in table lllll and 
figured Compared to the previous case of silver, only the 
values of the stiffness tensor were changed so most obser- 
vations previously made still apply. We recently reported 
on the experimental observation of the splitting of S2 for 
such nanoparticlesi^ There is an excellent agreement be- 
tween these measurements and the splitting calculated 
using the present approach. This strongly supports the 
validity of our approach and in particular the relevance of 
elastic anisotropy even in such very small nanoparticles. 

One notable difference between gold and silver con- 
cerns the case of the breathing mode. The isotropic 
breathing mode at 311.6 GHz belongs to the third Aig 
branch. Therefore it is tempting to assume that the 
anisotropic breathing mode lies on the same branch and is 
therefore mode i = 133 at 301.9 GHz. However, table IIIll 
reveals that mode i = 141 is a much better candidate 
for a breathing mode due to its larger volume variation 
and projection onto Sj. It turns out this is due to a 
strong mixing of the Aig branches as anisotropy is in- 
creased. Indeed, the isotropic modes S4, Se and Sg are 
split into various irreducible representations, one of them 
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u (GHz) 


i.r. 


7-8 


74.6 


Eg 


9-11 


76.9 


T2u 


12 


109.0 


Aag 


13-15 


111.2 


Tiu 


16-17 


114.1 


Eu 


18-20 


120.5 


T2g 


21-23 


140.2 


T2g 


24-26 


141.8 


Tig 


27-28 


154.0 


Eg 


29-31 


158.9 


T2u 


32-34 


168.2 


T2u 


35-37 


169.7 


Tiu 


43 


182.0 


-Lg 


118 


285.8 




133 


301.9 


Aig 


141 


310.1 


Aig 


183 


341.5 


Aig 



Squared Lamb projections 
0.996 Sa 
0.939 
0.996 Tl 
0.516 S3 
0.977 
0.837 S2 
0.773 Tl 
0.433 



0.824 S2 
0.698 Ti 
0.617 S3 
0.484 S; 



+ 0.002 
-I- 0.053 
+ 0.002 
+ 0.466 
+ 0.016 
+ 0.146 
+ 0.162 
+ 0.356 
+ 0.126 
-I- 0.259 
+ 0.186 
+ 0.456 



Si 
Si 

si 

s; 

n 

si + 
si + 
si + 
si + 

Ti -I- 

sl + 



0.987 Si + 0.005 Si + 
0.594 Sl + 0.218 Sg + 
0.866 Sl + 0.049 Sg -f 
0.806 Sl + 0.133 Sg + 
0.378 Sg + 0.319 sl + 



5V 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

0.1 

1.2 

0.3 

3.9 

1.4 



TABLE IIL Characteristics of the modes of a mono-domain 
sphere of Au with R = 5 ran. i is the mode index. The two 
largest projections onto the Lamb modes are shown. 



being Aig. Unlike the previous case of silver, the elastic 
constants of gold result in four Aig branches coming from 
Sj, S4, S| and S| being in the same frequency range. As 
a result, mode i = 141 which is the mode having the 
strongest projection onto Sj and also has a large volume 
variation does not lie on the same branch as Sj. Figure [3] 
shows the variation of the squared projections onto Sj of 
three of these Aig branches as the elastic anisotropy is 
varied and confirm the mixing discussed above. 

Similar but less pronounced mixings are observed in all 
cases. The density of modes (i.e. the number of modes 
per unit frequency) increases with frequency. As a result, 
the probability of mixings increases with frequency or i. 
This explains why the lowest frequency modes such as 
the lowest Eg ones are almost pure isotropic modes (83). 
However, the lowest T2g branches which are issued from 
the same S2 and from T| are already significantly mixed 
by anisotropy for both gold and silver nanoparticles. This 
is clearly evidenced in figure [3] too. 






isotropic 



anisotropic 



FIG. 2; (Color online) Evolution of the frequency of Raman 
active modes of a mono-domain gold sphere (radius 5 nm) 
with varying elastic anisotropy up to the frequency range of 
the breathing mode. Raman inactive modes are plotted with 
black lines. Eg branches with lines and crosses (red online), 
T2g branches with lines and triangles (red online) and Aig 
branches with lines and full circles (blue online). 



scattering by Sq and S2 modes is allowedii and has a sig- 
nificant Raman scattering cross-sectior^ in the isotropic 
case. Anisotropic nanoparticles obey the same rules and 
one might qualitatively estimate their Raman intensity 
by using the projections onto the same isotropic modes. 
As a result a significant scattering intensity is expected 
for the lowest Eg and T2g modes of both silver and gold 
while the second T2g mode should have a small Raman 
cross-section due to its small projection onto S2. The 
volume mechanism enables a measurable Raman inten- 
sity for the Aig modes having a significant projection 
onto So, i.e. for modes z = 75 and i = 141 for silver 
and gold respectively. The Raman intensity for other 
Aig modes should be very weak and probably not mea- 
surable in practice. 



2. Raman scattering efficiency for metallic cubic materials 

Low-frequency Raman scattering from silver and gold 
nanoparticles has attracted a lot of attention during the 
last few decades. The scattering mechanisms have been 
identified^^' and one might wonder how the anisotropic 
considerations detailed in this work fit into this picture. 
The inelastic light scattering process for such nanoparti- 
cles is mediated by the dipolar plasmon. As a result only 



3. Other cubic materials 

For cubic materials, the degree of elastic anisotropy 
is quantified by the Zener anisotropy ratio, A = 
2Cii/{Cii — C12). A=l for an isotropic material. Sil- 
ver and copper both have v4 ~ 3 and also had the largest 
mode splittings that we found. For comparison, we men- 
tion results for some less anisotropic materials. For a 
silicon nanosphere with radius 5 nm, the Lamb mode 
S2 splits into two frequencies at 386.6 GHz (Eg) and 
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FIG. 3: (Color online) Evolution of the squared projection of 
several modes of a mono-domain gold sphere (radius 5 nm) 
with varying elastic anisotropy. The bottom plot shows the 
squared projections onto S2 of the Eg branch coming from S2 
(line with circles, blue online) and the Tag branches coming 
from S2 (line with triangles, red online) and T3 (black line). 
The upper plot shows the squared projections onto Sj of the 
Aig branches coming from Sg (line and triangles, red online). 
So (black line) and S4 (line and circles, blue online). 



485.9 GHz (T2g). For a germanium nanosphcre, the same 
mode is split into 235.1 GHz (Eg) and 295.3 GHz (Tzg). 
Using the standard deviation as a simple measure of the 
frequency splitting A (and therefore neglecting the ef- 
fect of mixings), we obtain A =23, 22, 10 and 11% for 
Ag, Au, Si and Ge respectively. For all these commonly 
studied materials, this sphtting cannot be neglected. 



B. Spherical nanocrystals with tetragonal 
crystallinity 

Nanospheres with tetragonal crystallinity have a lower 
symmetry than the previous nanospheres with cubic crys- 
tallinity. The corresponding point group is D4h. As a 
result, more degeneracy lifting occurs. This results in 
a splitting of the infrared active Si mode into A2U and 
Eu and four branches starting from the S2 mode (Aig, 
Big, B2g and Eg). The relative order of these branches 
depends on the stiffness tensor. Two nanospheres made 
of Ti02 with radius 5 nm will be considered using the 
parameters from Ref. [2^. One of them has the anatase 
crystal structure (table lTV)) and the other the rutile struc- 
ture (t able fVl and figured]). We used these calculations for 
the anatase crystal structure in a recent worki to model 
the inelastic scattering of neutrons. Both have tetragonal 
symmetry but different elastic parameters. This results 
in different relative positions of the frequencies of the 
modes coming from a given isotropic mode. 
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1/ (GHz) 


i.r. 
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283.5 


Aig 


8-9 


284.0 


Eu 


10 


298.5 


Aiu 


11 


304.3 


Biu 


12-13 


316.1 


Eg 


14 


326.9 


Bag 


15 


330.9 


Bau 


16 


381.5 


Big 


17 


432.4 


Aau 


18 


443.3 


Bag 


19-20 


445.7 


Eu 


21 


447.6 


Aau 


22-23 


455.5 


Eg 


24-25 


463.4 


Eu 


26 


470.7 


Aag 


27 


473.9 


Big 


28 


476.9 


Biu 


29-30 


480.7 


Eg 


31 


510.9 


Bau 


32-33 


537.1 


Eu 


34 


573.8 


Aig 


35 


580.8 


Ai„ 
^g 


55 


678.4 


Aig 


67 


743.1 


Aig 


69 


751.6 


Aig 



Squared Lamb projections 
0.977 Sa 
0.931 T^ 
0.999 T^ 
0.998 T^ 
0.996 S^ 
0.992 Sa 
0.967 T^ 
0.965 Sa 
0.494 Si 
0.950 Tl 
0.612 Si 
0.537 S3 
0.957 T3 
0.758 S3 
0.989 Tl 
0.865 T3 
0.987 S3 
0.941 T3 
0.946 S3 
0.826 S3 
0.974 si 



0.910 Sl 



+ 0.016 
+ 0.054 
+ 0.000 
-I- 0.001 
+ 0.002 
+ 0.008 
+ 0.027 
+ 0.033 
+ 0.436 
+ 0.028 
+ 0.323 
+ 0.448 
+ 0.015 
+ 0.220 
-f- 0.005 
+ 0.060 
+ 0.008 
+ 0.031 
+ 0.026 
+ 0.095 
+ 0.009 
+ 0.030 



sl + 

Tl + 
Si + 
Tl + 

Tl + 
Si + 
Tl + 
Si + 
Si + 
Si + 

si + 
si + 
si + 

Tl + 
Si + 
Tl + 
Si + 
Tl + 
Si + 

si + 



0.885 Si + 0.082 Tj 
0.803 So -I- 0.156 Tl 
0.675 T^ + 0.162 So 



SV 



0.6 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.8 

0.2 

4.0 

1.8 



TABLE IV; Gharacteristics of the modes of a mono-domain 
anatase TiOa sphere with R = 5 nm are shown, i is the 
mode index. The two largest projections onto Lamb modes 
are shown. 



Using the same measure as before, the frequency split- 
ting for the S2 modes is A = 10 and 20% for the anatase 
and rutile structure respectively. 

Due to the lowering of the symmetry compared to the 
cubic case, the breathing mode is also susceptible to more 
mixings. For anatase Ti02, the Sj mode mixes mainly 
(mode i = 67) while for rutile Ti02 it mixes 



with Tg 



mainly with S2 (mode i — 41). 



C. Spherical nanocrystals with hexagonal 
crystallinity 

As a last example of the influence of elastic anisotropy, 
we focus on nanospheres made of crystals having hexago- 
nal symmetry, namely CdSe (table IVTl and figure E] using 
the elastic constants from Ref . 1271) . Co (table IVlTl using 
the elastic constants from Ref. (28|) and ZnO (elastic con- 
stants from Ref. 18). The associated point group is Dgh- 
The most notable difference compared to previous sym- 
metries is that some "accidental" degeneracy exists for 
such systems: for each Big (Biu) vibration there is a B2g 
(B2u) vibration having the same frequency. 

Using the same measure as before, the frequency split- 
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i 


u (GHz) 


i.r. 


7 


288.1 


Big 


8 


296.9 


B2u 


9-10 


355.8 


Eu 


11 


419.8 


Aiu 


12-13 


444.1 


Eg 


14 


453.1 


Aig 


15 


481.0 


Biu 


16 


499.4 


Big 


17-18 


531.7 


Eu 


19-20 


537.1 


Eg 


21 


539.4 


A2e 


22 


555.0 


Bag 


23 


558.9 


A2u 


24 


600.0 


B2u 


25-26 


622.3 


Eu 


38 


744.0 


Ala 


39-40 


744.2 


Eu 


41 


775.8 


Ai„ 


50 


836.1 


Aig 


62 


911.0 


Aig 


81 


1019.1 


Aig 



Squared Lamb projections 
0.978 Sa 
0.943 
0.693 
0.981 
0.927 
0.955 Sa 
0.965 
0.715 Tl 
0.402 Si 
0.813 Tl 
0.465 Tl 
0.997 Sa 
0.969 Sl 
0.845 S3 
0.575 S3 



+ 0.013 
+ 0.045 
+ 0.150 
+ 0.013 
+ 0.062 
+ 0.035 
4- 0.033 
-I- 0.210 
+ 0.288 
+ 0.071 
+ 0.327 
+ 0.001 
+ 0.016 
+ 0.095 
+ 0.247 



sl + 
sl + 

Ti + 
T^^ + 

50 + 

51 + 

Si + 

+ 
S^ + 
Si + 

si + 
sl + 

Ti + 
Sl + 



0.739 Si + 0.157 Sj + 
0.725 Sl + 0.178 Sl + 
0.631 Sl + 0.269 Si + 

0.439 si -f- 0.351 Si + 

0.739 Si + 0.169 Si + 

0.649 + 0.110 Sl + 



SV 



0.0 
0.0 
0.0 
0.0 
0.0 
0.8 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

1.4 
0.0 
4.1 

1.6 

0.9 

0.5 



TABLE V: Characteristics of the modes of a mono-domain 
rutile TiOa sphere with R = 5 nm are shown, i is the mode in- 
dex. The two largest projections onto Lamb modes are shown. 
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FIG. 4; (Color onhne) Evolution of the frequency of Raman 
active modes of a mono-domain rutile Ti02 sphere (radius 
5 nm) with varying elastic anisotropy up to the frequency 
range of the breathing mode. Raman inactive modes are plot- 
ted with black lines, Eg branches with lines and crosses (red 
online), Big branches with lines and asterisks (red online), 
B2g branches with lines and triangles (red online) and Aig 
branches with lines and full circles (blue online). 



ting for the S2 modes is A =9, 11 and 5% for the CdSe, 
Co and ZnO respectively. ZnO is therefore the most elas- 
tically isotropic amongst the materials studied in this 
work. 

Regarding the breathing mode, for CdSe it is hardly 
mixed with other modes due to anisotropy and results in 
mode i = 77 which has a strong projection onto Sj and 
a large volume variation. For the nanosphere made of 
cobalt, there is a strong mixing with the and Si modes 
and the anisotropic mode with the largest projection onto 
Sq and the largest volume variation is mode z = 49 which 
is on the branch coming from S4 at 466.7 GHz and not 
on the branch starting from Sq at 536.8 GHz. 

The infrared active S} mode which was recently 
observed^ for CdSe nanoparticles is split by the lower- 
ing of symmetry into and Eiu modes at 181.9 and 
181.3 GHz respectively. This frequency splitting is very 
small but these modes mix with the branches having 
identical irreducible representations coming from S|. Ne- 
glecting this mixing, the isotropic is in excellent agree- 
ment with the anisotropic description. 



D. Non-spherical nanocrystals 

In order to illustrate the usefulness of the same nu- 
merical tools for a different source of anisotropy, we con- 



sider in the following the lowering of the symmetry due 
to the shape of the nanoparticles. We start with a mi- 
nor change of the shape as the nanosphere is transformed 
into a spheroid having an isotropic elasticity. Then we 
consider faceted nanoparticles with elastic anisotropy. 



1. Spheroids with isotropic elasticity 

Let us consider a spheroid made of silver with a degen- 
erate semi-axis R = 5 nm and a varying non-degenerate 
semi-axis Rz- We assume the elasticity to be isotropic. 
The point group associated with such a spheroid is there- 
fore Dooh- The frequencies, irreducible representations 
and volume variations for a spheroid with Rz — 10 nm 
are presented in table IVIIII The branches obtained for 
varying Rz are presented in figure [51 

Let us now focus on the effect of the spheroidal defor- 
mation on the S2 modes. Looking at the displacements 
corresponding to these modes for small deviations from 
the sphere, it is possible to understand the frequency vari- 
ations. The lowest Aig mode corresponds to a stretch- 
ing along the z direction accompanied by a shrinking in 
the xy plane. Therefore it can be seen as a vibration 
confined along the z direction and its frequency varies 
roughly as 1/Rz- The E2g vibrations correspond to a 
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7 

8-9 
10-11 
12-13 
14-15 

16 
17-18 

19 

20 
21-22 
23-24 
25-26 
27-28 
29-30 
31-32 

33 

55 
56 
57 

77 



V (GHz) 



119.4 
123.2 
126.8 
132.6 
132.8 
158.6 
181.3 
181.9 
186.1 
189.2 
192.1 
193.5 
198.1 
208.0 
214.7 
226.7 

273.8 
284.9 
296.0 

320.4 



Aiu 

E2u 

Eig 

E2g 

Eiu 
Aig 
Eiu 

A2u 

Big-fB2g 

E2u 

Eig 

Blu+B2u 
E2g 

Eiu 

A2u 

Aig 

A2g 

Aig 
Aig 



Squared 
1.000 
0.998 
0.999 Sa 
0.998 Sa 
0.971 
0.999 Sa 
0.894 Si 
0.789 S} 
0.999 
0.994 % 
0.998 S3 
0.980 T\ 
0.996 S3 
0.947 T3 
0.870 Sj 
0.789 S3 



Lamb projections 
+ 0.000 + 
+ 0.001 Sj + 
+ 0.001 + 
+ 0.002 + 
+ 0.023 + 
+ 0.001 si + 
+ 0.103 Si^ 4- 
-)- 0.206 S3 + 
0.001 T} + 
+ 0.004 Si + 
+ 0.001 + 
+ 0.014 Si + 
+ 0.004 Ti + 
+ 0.042 si 
-h 0.099 si -I- 
-f 0.206 Sl -\- 



0.926 Si + 0.060 Si -I- 
0.999 T\ 0.001 Ti + 
0.925 Si + 0.059 Si + 

0.985 So + 0.011 Si 4- 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

0.3 
0.0 
0.2 

4.5 



350 



TABLE VL Characteristics of the modes of a mono-domain 
wurtzite CdSe sphere with R = 5 nm are shown, i is the 
mode index. The two largest projections onto Lamb modes 
are shown. 



i 


u (GHz) 


i.r. 


7-8 


229.8 


Eau 


9 


232.2 


Am 


10-11 


240.7 


Eag 


12-13 


246.6 


Eig 


14-15 


251.4 


Eiu 


16 


316.1 


Aig 


17-18 


334.7 


Eiu 


19 


346.4 


Aau 


49 


497.8 


Aig 



Squared Lamb projections 

T 



T^ + 0.001 Ti + 
1.000 T^ -I- 
0.997 Sa + 
1.000 Sa 4- 
0.942 Ti + 
0.996 Si + 
0.919 Si + 
0.828 si -f- 



0.000 Ti 
0.003 Ti 
0.000 Ti 
0.029 Si 

0.003 si 

0.063 Si 
0.167 Si 



0.958 Si + 0.042 Si 



5V 



0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

4.5 



TABLE VH: Characteristics of the modes of a mono-domain 
Co sphere with R — 5 nm are shown, i is the mode index. 
The two largest projections onto Lamb modes are shown. 



stretching in the xy plane without changes along the z 
axis and therefore their frequencies hardly changes with 
Rz- The Eig vibrations correspond to a stretching in the 
xz and yz plane without changes along the z axis and 
therefore their frequencies vary with a slower eccentric- 
ity dependence than the previous mode. These rough 
approximations are in agreement with the dependence 
observed in figure [S] 

Using perturbation theory^^^ it is possible to obtain 
more accurate expressions for the frequencies of these 
three branches. The exact variations for \Rz — R\ R 
are 17(1 + 4/3/21), ^{1 - 2(3/21) and ^{1 - 4/3/21) for the 
Eig, E2g and Aig modes respectively where (i = 2{Rz — 






isotropic 



anisotropic 



FIG. 5; (Color online) Evolution of the frequency of Raman 
active modes of a mono-domain wurtzite CdSe sphere (radius 
5 nm) with varying elastic anisotropy up to the frequency 
range of the breathing mode. Raman inactive modes are plot- 
ted with black lines, Eig branches with lines and crosses (red 
online) , Eag branches with lines and triangles (red online) and 
Aig branches with lines and full circles (blue online). 



R)/{Rz + R) and il is the frequency of modes Si for a 
spherical particle having the same volume. Note that 
the length of the degenerate semi-axis is constant in this 
work and therefore the volume varies linearly with Rz- 
Using fl = wgi y^R/Rz results in expressions which are 
in very good agreement with figure [6] close to Rz = R- 

Regarding the "breathing" mode, the picture gets very 
complicated when Rz/ R differs significantly from 1. Sj 
is on the 4**^ Aig branch. Therefore the mode of the 
spheroid with Rz/R = 2 on the same branch is mode 
j = 64. However, due to the mixing with neighboring 
Aig branches, the mode with i = 90 is a better candi- 
date since its volume variation is much larger. Thanks 
to this additional property of the symmetric modes, it 
is therefore possible to follow the breathing mode. How- 
ever, because it is not possible to project onto the Lamb 
modes of an isotropic sphere due to the different shapes, 
only the branches and the anti-crossing patterns between 
branches having the same irreducible representation can 
help tracking qualitatively the other S and T modes. 

These calculations are targetted at interpreting ex- 
perimental results on silver nanoparticles such as 
nanocolumns."^*^ Compared to previous calculations us- 
ing the FEMS method presented before, '^"'^ the current 
approach enables a more complete description of the dif- 
ferent vibrations besides being faster and more acurate. 
In particular, the variation of volume is very efficient 



9 



i 
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SV 


7-8 


62.4 


Eiu 


0.0 


9 


71.1 


A2u 


0.0 
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90.8 
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13 
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0.0 


14-15 


127.0 


Eig 


0.0 


16-17 


136.3 


J_/2g 


0.0 


34 


185.2 


i-iig 


0.2 


57 


244.6 


Ai 


1.0 


64 


253.3 


Aig 


0.0 


90 


294.8 


Aig 


4.5 


110 


318.8 


Aig 


0.5 



TABLE VIII: Characteristics of modes of a spheroid of dimen- 
sions R = 5 nm, Rz = 10 nm made from elastically isotropic 
{i.e. not mono-domain) Ag. i is the mode index. 




0.5 1 1.5 2 

R^/R 



FIG. 6: Reduced frequency of the vibrations of an elastically 
isotropic silver spheroid as a function of its aspect ratio. R 
is the length of the degenerate semi-axis and Rz is the length 
of the non-degenerate one. a^gi is the frequency of the S2 
mode of the sphere of radius R. Raman inactive modes are 
plotted with black lines, Eig branches with lines and crosses 
(red online), E2g branches with lines and triangles (red online) 
and Aig branches with lines and full circles (blue online). 



in showing which modes should be observed by time- 
resolved pump-probe experiments.-^ It is interesting to 
note that some work o^'^'^ concern aligned nanocolumns 
which results in interesting depolarization rules for the 
Raman peaks. For non-aligned nanoparticles, these rules 
are the same as those used routinely for an ensemble 
of molecules, i.e. all the Raman active modes produce 
completely depolarized Raman peaks except for the Aig 
vibrations which have a "polarized" scattering for which 
the Raman peak is more intense when the polarizations 



of the incident and scattered photons are parallel. For 
oriented nanoparticles, this rule doesn't hold and the an- 
gles of the incident and scattered photons with respect 
to the axis of symmetry of the nanoparticles have to be 
taken into account. 

Recent time-resolved pump-probe femtosecond exper- 
iments for single gold nanoparticles'^^ also demonstrate 
the need for such a model. In particular, the observa- 
tion of a peak close to the S2 frequency is reported for 
spheroids but not for spheres. This is in agreement with 
the features reported here for silver, namely that the Aig 
mode coming from the S2 mode has a non-zero volume 
variation enabling it to be observed in such an experiment 
unlike the S2 mode of a sphere. It is also worth noting 
that the shape of the peak attributed to "the breathing 
mode" in this work for dumbbells looks quite compli- 
cated. Such nano-objects have the same symmetry as 
spheroids. As a result, in most cases there is no such 
thing as a "breathing mode" but rather a set of Aig vi- 
brations in a relatively narrow frequency range having a 
significant volume variation. This is due to the fact that 
the Aig branch coming from the Sj vibration mixes with 
all the other Aig branches which come from all the Sg 
vibrations with even £. Modeling such gold dumbbells is 
beyond the scope of this paper, but doing so would en- 
able a more detailed understanding of the experimental 
results, especially for such single particle measurements 
for which the external shape of the nanoparticles can 
be obtained from SEM images. The relatively large size 
of these nanoparticles prevents them from being single- 
domain and justifies the use of the isotropic elastic ap- 
proximation. 



2. anisotropic gold polyhedra 

Observation of faceted nanoparticles using electron mi- 
croscopy is quite common. Such facets can be thought of 
as a signature of the inner crystal structure and therefore 
as an indication of the elastic anisotropy.'^^ To quantify 
the importance of the shape on the vibrations, we cal- 
culated the frequencies and irreducible representations 
for some polyhedra and the results are presented in ta- 
ble |TXl The crystal lattice is oriented with respect to the 
shape so that the [100] planes correspond to the square 
faces for the cuboctahedron and to the octagonal faces 
for the truncated cuboctahedron. There is no lowering 
of symmetry associated with these shapes compared to 
the case of a mono-domain spherical gold nanoparticle. 
The relevant point group is then . While table IIXI 
clearly shows that the lowest frequencies change with the 
shape, the frequencies of the lowest Eg and T2g modes 
are hardly affected. Since the degeneracies of the So 
and Si modes are not lifted, all the modes which are 
observable by Raman scattering, infrared absorption of 
time-resolved pump-probe experiments are not sensitive 
to these changes of shape. 
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i.r. 


sphere 


cuboctahedron 


truncated cuboctahedron 


Eg 


74.6 


74.5 


74.5 


T2u 


76.9 


73.4 


75.5 


Aag 


109.0 


81.6 


91.7 


Tiu 


111.2 


108.7 


110.7 


Eu 


114.1 


102.3 


107.9 


T2g 


120.5 


121.3 


121.7 


T2g 


140.2 


128.0 


133.8 



TABLE IX: Frequencies and irreducible representations of 
the lowest frequency vibrations of elastically anisotropic gold 
nanocrystals having different shapes. The volume of the dif- 
ferent nanocrystals is the same as that of a sphere of radius 
5 nm. The frequencies of the Raman active modes (Eg and 
T2g) are almost unaffected by the shape. See the text for the 
orientation of the crystal lattice with respect to the polyedra. 

IV. MULTIPLE-DOMAIN NANOPARTICLES 

Several attempts have been made in the past either to 
fit low frequency Raman spectra or to determine the size 
distribution of the nanoparticles inside a sample using 
the shape of the low-frequency Raman peak. Both ap- 
proaches always rely on the validity of the isotropic model 
by Lamb and on the predominance of the size distribu- 
tion, the coupling with a surrounding matrix and the 
electron-vibration coupling to fit the broadening of the 
peaks. However the distribution of internal structures of 
multiple-domain nanoparticles also results in inhomoge- 
neous broadening. 

It is in principle possible to model the vibrations 
of a multiple-domain nanoparticle using the numerical 
method of Visscher et at. However this requires a com- 
plete description of the position of the domain boundaries 
and the orientations of the crystal lattice. For ensemble 
measurements with the nanoparticles having a variety of 
different internal structure, a lot of calculations would be 
required. Otherwise, using such an approach for a sin- 
gle internal structure is justified only if all the studied 
nanoparticles are identical for ensemble measurements 
or if the inner structure of a nanoparticle studied in a 
single particle measurement is perfectly known. Due to 
these latter two conditions having never been met un- 
til now and also to the additional complexity of mod- 
elling a multiple-domain nanoparticle, we suggest using 
the isotropic approximation to describe an ensemble of 
multiple-domain nanoparticles. What this means is that 
no nanoparticle behaves exactly as an isotropic nanopar- 
ticle, but the isotropic approximation gives an average 
value due to the nanoparticles having essentially ran- 
dom domain structures. Ensemble measurements should 
therefore show features associated with these average fre- 
quencies with some inhomogeneous width due to the in- 
homogeneous distribution of domains inside the popula- 
tion of nanoparticles. 

We calculated the vibrations of an icosahcdron made 



of gold and silver assuming the isotropic approximation 
to be valid in that case. The corresponding point group 
for such a nanoparticle is Ih. It is interesting to note 
that in this case there is no degeneracy lifting for the 
modes So, Si and S2 whose irreducible representation in 
the new system is Ag, Tiu and Hg respectively. As a 
result, there is hardly any frequency difference compared 
to the case of a sphere having the same volume. A real 
icosahedron made of an anisotropic material but having 
the same symmetry due to the presence of twins would 
have no degeneracy lifting either for the same modes and 
we expect almost the same frequencies too. 



V. CONCLUSION 

We have presented mode frequencies and irreducible 
representations for homogeneous continuum nanoparti- 
cles using a standard numerical method which can handle 
arbitrary shape and anisotropic elasticity. The classifica- 
tion by irreducible representation makes it possible to 
label many modes as either Raman or infrared inactive. 
We have been able to go beyond this to provide some 
tools to make qualitive estimates of the Raman inten- 
sity of potentially Raman active modes. These tools are 
well-suited for the interpretation of experimental results 
obtained with vibrational spectroscopiesi^^ 

This approach fills a large gap in current works since 
up to now it was necessary to choose between a sim- 
plified spherical isotropic model where the underlying 
physics was simple enough but where the accuracy was 
challenged by recent experimental results or a numerical 
approach which provides accurate frequencies but which 
is inefficient in practice due to the lack of tools to dis- 
tinguish the relevant vibrations without complex simula- 
tions of spectra. 

It is clear that additional theoretical work is required 
in order to quantitatively predict Raman intensities for 
elastically anisotropic modes of spheres. However, the la- 
beling of modes via their Lamb mode parentage by con- 
tinuous variation of the elastic isotropy and projections 
for spherical nanoparticles, is a powerful descriptive and 
semi-quantitative tool for understanding what is going on 
with elastic anisotropy. The significant frequency split- 
tings obtained in this work call into question the validity 
of the isotropic approximation for the case of multiple- 
domain nanoparticles. Such systems are very relevant 
experimentally but their vibrations remain largely unad- 
dressed. 
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